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Abstract 

The high-pressure structural behavior of lanthanum monochalcogenides is investigated by theory 
and experiment. Theory comprises density functional calculations of LaS, LaSe and LaTe with the 
general gradient approximation for exchange and correlation effects, as implemented within the 
full-potential linear muffin-tin orbital method. The experimental studies consist of high-pressure 
angle dispersive x-ray diffraction investigations of LaS and LaSe up to a maximum pressure of 41 
GPa. A structural phase transition from the NaCl type to CsCl type crystal structure is found 
to occur in all cases. The experimental transition pressures are 27-28 GPa and 19 GPa, for LaS 
and LaSe, respectively, while the calculated transition pressures are 29 GPa, 21 GPa and 10 GPa 
for LaS, LaSe and LaTe, respectively. The calculated ground state properties such as equilibrium 
lattice constant, bulk modulus and its pressure derivative, and Debye temperatures are in good 
agreement with experimental results. Elastic constants are predicted from the calculations. 
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I. INTRODUCTION 



The lanthanum monochalcogenides belong to the wide class of binary rare-earth 
monochalcogenides with the NaCl-type structure, which has been intensively studied be- 
cause of their interesting physical properties including complex magnetic structures, pres- 
sure induced insulator-metal transitions, anomalous valence fluctuations, and unusual Fermi 
surface properties.-^ The trivalent lanthanum monochalcogenides are superconductors near 
1 K, and the superconducting transition temperature as well as the electronic specific heat 
coefficient increases from the monosulfide to the monotelluride, whereas the Debye temper- 
ature decreases from LaS to LaTe^. 

In this work we explore the high pressure behavior of the La chalcogenides by experiment 
and theory. High-pressure x-ray diffraction experiments are conducted on LaS and LaSe 
to obtain equations of states including discontinuous structural transitions from the NaCl 
structure (Space group Fm3m, Z=l, also called the Bl structure in the following) to the 
CsCl structure (Space group Pm3m, Z=4, also termed the B2 structure in the following). 
Density functional calculations are performed to compare basic theoretical predictions with 
the measurements. The present experimental study of LaS extends our previous study on 
U x Lai_ x S compounds^. For LaSe no high-pressure results have been reported previously, 
while high-pressure experiments on LaTe report a structural transition from Bl to B2 at 
around 7 GPa£. 

A second objective of the present work is to investigate the elastic properties of the 
lanthanum chalcogenides, for which there are no experimental results available. Elastic 
constants are derived from total energy variations with applied strains. We compare our 
results with recent theoretical work& and find good agreement. 

The systematics of the electronic band structure through the lanthanum and cerium 
chalcogenide series was studied with angle-resolved photoemission spectroscopy^. The 
Fermi surface properties of lanthanum monochalcogenides have been studied by de Haas- van 
Alphen effect measurement s*^. From the theoretical side the energy band structure, su- 
perconductivity, surface electronic structure, optical and magneto-optical spectra as well 
as structural stability of lanthanum monochalcogenides have been investigated by sev- 
eral authors using the local spin density approximation (LSDA) as well as the LSDA+U 
approximation^^^^ 1 ^ 1 ^ 1 ^. A few papers also focused on the calculation of second and 



third order elastic constants using the short range repulsive potential methocU^^. A 
few experimental studies such as point contact spectroscopy—, reflectivity^ and phonon 
spectra^ have been reported for LaS. The magneto-optical properties of lanthanum chalco- 
genides have been investigated experimentally,— 1 ^ 1 ^ where a non-zero Kerr effect could be 
observed in these Pauli paramagnets by applying an external magnetic field. 

The rest of the paper is organized as follows. In section 2 we discuss the computational 
and experimental details of our work. The calculated ground state properties, elastic con- 
stants as well as the experimental and theoretical results for the high-pressure behavior and 
structural transitions are presented in Section 3, and the conclusions are given in Section 4. 



II. COMPUTATIONAL AND EXPERIMENTAL DETAILS 



A. The electronic structure method 



In this work we have used the all-electron full-potential linear muffin tin orbital (FP- 
LMTO) method^ to calculate the total energies and the basic ground state properties. Here 
the crystal is divided into two regions: non-overlapping muffin-tin spheres surrounding each 
atom and the interstitial region between the spheres. We used a double k spdf LMTO basis 
(each radial function within the spheres is matched to a Hankel function in the interstitial 
region) for describing the valence bands. The following basis set was used in the calculations: 
La(5s,6s,5p,5d,4f), S(3s,3p,3d), Se(4s,4p,3d,4d), Te(5s,5p,4d,5d). The exchange correlation 
potential was calculated within the generalized gradient approximation (GGA) scheme^. 
The charge density and potential inside the muffin-tin spheres were expanded in terms of 
spherical harmonics up to l max =Q, while in the interstitial region, they were expanded in 
plane waves, with 6566 (energy up to 109.80 Ry) waves being included in the calculation. 
Total energies were calculated as a function of volume, for a (18 18 18) k-mesh, corresponding 
to 195 k- vectors in the irreducible wedge of the Brillouin zone, and the results fitted to the 
Birch equation of stated to obtain the ground state properties. 

The elastic constants were obtained from the variation of the total energy under volume- 



conserving strains, as outlined in Ref. 
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B. Experimental details 



Experiments were carried out at the European Synchrotron Radiation Facility (ESRF), 
on the ID30 beamline dedicated to high pressure diffraction experiments, in the angular 
dispersive mode, using a double focused monochromatic beam at A=0.3738 A. Loading of 
LaS and LaSe powders were performed in Le Toullec-type diamond anvil cells (DAC), using 
nitrogen, argon or silicone oil as the pressure transmitting media. The pressure inside the 
cell was determined via the ruby fluorescence method^ from a small ruby ball mounted 
alongside the LaX sample. The sample to detector distance was calibrated before each 
set of experiments by means of a standard silicon powder sample. Diffraction images were 
captured with a FASTSCAN— image plate detector, and processed using the ESRF Fit2D 
program^ to provide data in a format suitable for Rietveld analysis using the FULLPROF 
program^. The lattice constants at ambient pressure have been determined as ao=5.852 and 
a =6.067 A for LaS and LaSe, respectively, in good agreement with previous experimental 
data given by Ref. Il6 



III. STRUCTURAL, ELASTIC AND HIGH PRESSURE STUDIES 
A. Ground state and Elastic properties 

The calculated ground state properties such as equilibrium lattice constant, bulk modulus 
and its pressure derivative are given in Table 1. The calculated equilibrium lattice constants 
are generally overestimated by ~ 0.5 % compared to experiment, which is an improvement 
compared to earlier calculations^ 1 ^. Similarly, for the bulk modulus, the agreement between 
theory and experiment has improvecr^£. It is interesting to note that the bulk moduli of the 
lanthanum chalcogenides are quite similar to those of the neighboring cerium chalcogenides, 
which are characterized by a localized / stated B = 82 GPa, 76 GPa and 58 GPa, for CeS, 
CeSe and CeTe, respectively. Similarly, the calculated bulk moduli of the praseodymium 
chalcogenides are B = 89 GPa, 78 GPa and 57 GPa, for PrS, PrSe and PrTe, respectively, 
when Pr is represented in a trivalent configuration^. These similarities of numbers corrob- 
orate the practice of taking the lanthanum chalcogenides as non-magnetic reference systems 
for the later rare-earth chalcogenides. 

The calculated elastic constants of LaS, LaSe and LaTe are listed in Table 2, where we also 
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compare to the recent LAPW calculation of Ref. |6|. The two theoretical approaches find very 
similar Cu value for LaS and LaSe, while for LaTe this parameter is ~ 13% smaller in the 
present calculations than found by the LAPW method. The present Cyi and C44 parameters 
are in all cases lower than found in Ref. |(3, for LaSe C12 with almost a factor of 2 difference. 
Our calculated shear moduli G are 19% to 40% larger than the ones calculated with LAPW. 
For Young's modulus E, our values are 16% to 35% higher, and finally, for Poisson's ratio v 
our values are 8% to 17% lower than the LAPW results. The largest discrepancies between 
the two calculations occur for LaTe in all cases. Furthermore, we observe the general trend 
that the heavier the chalcogenide is, the softer is the compound. This trend is expected and 
easy to understand, since with increased chalcogen size, the lattice parameter increases and 
the valence orbitals become increasingly delocalized due to a higher number of nodes.— 

As previously stated, there is at present no experimental information regarding the elastic 
constants available for the lanthanum chalcogenides. The elastic constants of LaS are of 
roughly the same magnitude as those measured for the isoelectronic YS compound (Table 
2). 

Table 3 presents sound velocities and Debye temperatures, as derived from the calculated 
elastic constants^. The Debye temperatures have been determined experimentally from 
the low-temperature specific heat, and the agreement with our calculations is excellent (the 
calculated Debye temperatures are 0.5% to 6% higher than experiment), which can be taken 
as an indirect check on our calculated elastic constant values. 

B. High-pressure structural transitions 

The high pressure structural behavior of the lanthanum chalcogenides was studied both 
experimentally and theoretically. The calculated total energies as functions of relative vol- 
ume for LaS, LaSe and LaTe are shown in figure dj while experimental and theoretical 
pressure-volume relations are presented in figures 2, 4 and 5 for LaS, LaSe and LaTe, re- 
spectively. Figure 3 shows the diffraction spectra recorded for LaS and LaSe. It appears 
from the theory presented in figure 1 that with compression the B2 phase becomes more 
and more favorable, and eventually a transition from the Bl structure to the B2 structure 
occurs. From the common tangent of the Bl and B2 total energies the transition pressure is 
determined. The calculated transition pressure for LaS is 29.3 GPa with a volume collapse 



of 10.3%, which is in good agreement with the experimental transition pressure of 27-28 
GPa with a volume collapse of around 9.5% as shown in figure [2j A similar transition was 
predicted to occur in CeS around 24.3 GPa, in which the Ce ions are in a tetravalent stated 
In the case of PrS a similar transition is predicted to occur around 22 GPa, in this case with 
the Pr ions remaining trivalent^. The transition pressure of LaS is similar to that of the 
pnictogen group neighbour LaP, in which a transition from the Bl structure to a distorted 
B2 structure is observed experimentally around 24 GPa4^. 

Like LaS, LaSe also undergoes a Bl — > B2 transition, experimentally seen at a pressure 
of 19 GPa with a volume reduction of 10%, cf. figure HI for which the present theory finds 
concordant values of P t = 21 GPa with a volume collapse of 10.4%. In the case of CeSe a 
similar transition from Bl type to B2 type occurs around 20 GPa4^, in this case with the 
Ce ion remaining trivalent on both sides of the transition according to theory^. Even for 
PrSe a Bl — > B2 transition is predicted to occur around 12 GPa, wherein the Pr ions remain 
trivalent^. The transition pressure of LaSe is similar to that of the pnictogen neighbour 
LaAs, in which a transition from Bl structure to a distorted B2 type occurs around 20 
GPa^. 

As far as LaTe is concerned a transition from Bl type to B2 type is reported around 
7 GPa£, which agrees quite well with the calculated transition pressure of 9.7 GPa with a 
volume reduction of 10.4%, cf. figure 5. Both CeTe and PrTe undergo a transition from Bl 
type to B2 around 8±1 GPa4^ and 9±1 GPa^. When comparing the pnictogen neighbour 
LaSb with that of LaTe the pnitogen undergoes a transition from Bl type to a distorted B2 
type structure around 11 GPa— which is slightly higher than that of LaTe. 



IV. CONCLUSIONS 



By means of a combined theoretical and experimental study the high-pressure structural 
behavior of LaX (X=S, Se, Te) compounds have been investigated. Unlike the lanthanum 
monopnictides which show a transition from Bl to a distorted B2 structure^ 1 ^, the lan- 
thanum monochalcogenides exhibit a simple Bl — > B2 structural phase transition, similar 
to the one found for most of the lanthanide and actinide monochalcogenides studied up to 

For all three systems studied here, the volume collapse is around 10% (from both experi- 



ment and theory) at the Bl — > B2 transition. Let us compare this magnitude of the volume 
collapse to a very simple model in which we assume that at the transition point, the La-X 
bond length is restored to the ambient value. We expect such a model to underestimate the 
volume collapse, since the bond length at high pressure most probably will be a bit smaller 
than at ambient pressure. Our simple model gives the following relation between the relative 
volume V t B1 /V B1 just before the transition and the volume collapse 1 — V t m /V t B1 : 

yB2 a yBl 

1-1^ = 1- 4 v ° (1) 

For LaS, the experimental relative volume at the transition is 0.82, which would give a 
volume collapse of 6% according to this model. For LaSe, the corresponding numbers are 0.84 
and 8%. Thus, the volume collapse observed experimentally (and from density functional 
calculations) can be viewed as a partial restoration of the original bond length between 
the La and the chalcogenide. Based on the above discussion we conclude that the volume 
collapse observed in the lanthanum chalcogenides is entirely consistent with a simple picture 
of the transition in which the volume collapse is a consequence of the rearrangement of the 
atoms to a more close-packed structure while the total number of valence electrons remains 
unchanged. We have also checked the partial occupation number of spdf states across the 
structural transition and there is no appreciable change in the occupation numbers across 
the transition. This is in contrast to the situation in many Ce systems, in which the volume 
collapse is accompanied by derealization of an / electron.— 
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TABLE I: Calculated lattice constants in A, bulk modulus in GPa, and its pressure derivative B' 
of the lanthanum monochalcogenides in the Bl structure 



Compound 



Lattice Constant 



Bulk Modulus 



B' 



LaS 



LaSe 



LaTe 



Expt. this work 
Theory, this work 

Other theory 
Expt. this work 
Theory, this work 
Other theory 
Expt. 
Theory, this work 



5.852 
5.873 

5.812 a ,5.773 6 ,5.727 c ,5.895 e 
6.067 
6.091 
5.957 c ,6.126 e 
6.435 ft 
6.470 
6.255 c ,6.512 e 



89(3),83.6 d 
87.8 

97.74 a ,107 c ,81.53 e 
74(2) 
74.8 
97.74 c ,68.40 e 
60.6±2^,55 9 
59.4 



6.5(4) 
3.95 
4.67 e 
4.7(3) 
4.12 
4.28 e 

4.12 
4.96 e 



































a Ref 


12 


; 6 Ref 


14 


; c Rcf 
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; d Ref 
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; e Ref 
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; ^Ref 
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; b Ref. 


3; ft Ref. 
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TABLE II: Calculated elastic constants, shear modulus G, Young's modulus E - all given in GPa 
- and Poisson's ratio v for lanthanum monochalcogenides in the Bl structure at the theoretical 
equilibrium volume. For comparison, results of LAPW calculations (Ref. [6|) and experimental 



results for YS (Ref. 



40j) are included. 



Compound 


C u 


Cl2 


C44 


G 


E 


V 




LaS 


227.9 


18.0 


22.2 


55.3 


137.2 


0.240 


Present 




234 


23 


25 


46.6 


117.8 


0.26 


LAPW 


YS 


250 


20 


30 








Expt. 


LaSe 


201.6 


11.4 


15.7 


47.5 


117.5 


0.238 


Present 




203 


21 


22 


40.6 


102.5 


0.26 


LAPW 


LaTe 


158.7 


9.7 


7.9 


34.5 


86.8 


0.256 


Present 




171 


12 


8 


24.6 


64.1 


0.30 


LAPW 



1 1 



TABLE III: Calculated longitudinal, shear, and average wave velocity (vi, v s , and v m , respec- 
tively) in m/s, and the Debye temperature, 9d, in kelvin from the average elastic wave velocity for 
lanthanum monochalcogenides in the Bl structure at the theoretical equilibrium volume. 



Compound v\ v s v m 0£> 

LaS Present 5249 3072 3406 277.2 

Expt. a - - - 276 

LaSe Present 4626 2712 3006 233.6 

Expt. a - - - 231 

LaTe Present 3972 2273 2524 185.7 

Expt. a - ^ - - 175 

"Ref'l 



TABLE IV: Calculated and experimental transition pressures, Pt in GPa, and volume changes, 
given in %, for the Bl — ► B2 structural phase transition of lanthanum monochalcogenides. 



Compound 




Pi 




Volume Collapse 




Theory 




Expt. 


Theory Expt. 


LaS 


29.3 




27-28 


10.3 9.5 


LaSe 


21 




19 


10.4 10.5 


LaTe 


9.7 




r 


10.4 



RefH 



19 
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FIG. 1: (Color online) Calculated total energy vs. relative volume for (a): LaS, (b): LaSe, and (c): 



LaTe in the Bl and B2 structures. The common tangents mark the structural Bl — ► B2 transition. 
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FIG. 2: (Color online) Pressure as a function of relative volume for LaS in the Bl and B2 structures. 
Triangles: theory (GGA), circles: experiment. The experimental volume collapse at the Bl — > B2 
transition is marked. 
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. 3: X-ray diffraction spectra of (a): LaS and (b): LaSe in the Bl structure at P = GPa, 
5.8520 A for LaS, a = 6.0670 A for LaSe, and in the B2 structure right after the transition, at 
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FIG. 4: (Color online) Pressure as a function of relative volume for LaSe in the Bl and B2 
structures. Triangles: theory (GGA), circles: experiment. The experimental volume collapse at 
the Bl — ► B2 transition is marked. 
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FIG. 5: (Color online) Theoretical pressure as a function of relative volume for LaTe in the Bl and 
B2 structures. 



